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Abstract 



We consider the effects on big bang nucieosynthesis (BBN) of the radia- 
tive decay of a long-lived massive particle. If high-energy photons are emitted 
after the BBN epoch (t ~ 1 — 10 3 sec), they may change the abundances of 
the light elements through photodissociation processes, which may result in a 
significant discrepancy between standard BBN and observation. Taking into 
account recent observational and theoretical developments in this field, we re- 
vise our previous study constraining the abundance of the radiatively-decaying 
particles. In particular, on the theoretical side, it was recently claimed that 
the non-thermal production of 6 Li, which is caused by the photodissociation 
of 4 He, most severely constrains the abundance of the radiatively-decaying 
particle. We will see, however, it is premature to emphasize the importance 
of the non-thermal production of 6 Li because (i) the theoretical computation 
of the 6 Li abundance has large uncertainty due to the lack of the precise 
understanding of the 6 Li production cross section, and (ii) the observational 
data of 6 Li abundance has large errors. 
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I. INTRODUCTION 



Big-bang nucleosynthesis (BBN) is one of the most important tools to probe the early 
universe. Because it is very sensitive to the condition of the universe from 10 -2 sec to 
10 12 sec, we can indirectly check the history of the universe and impose constraints on a 
hypothetical particles by the observational light element abundances. 

There are a lot of models of the modern particle physics beyond the standard model, e.g. 
supergravity or superstring theory, which predict unstable massive particles with masses of 
(9(100 GeV), such as gravitino, Polonyi field, moduli, and so on. They have long lifetimes 
because their interactions are suppressed by inverse powers of the gravitational scale. Con- 
sequently, these exotic particles may decay at about the BBN epoch (T < lMeV). If the 
massive particles radiatively decays, the emitted high energy photons induce the electromag- 
netic cascade process. If the decay occurs after the BBN starts, the light elements would 
be destroyed by the cascade photons and their abundances would be changed significantly. 
Comparing the theoretically predicted light element abundances with the observational ones, 
we can impose constrains on the energy density, the mass, and the lifetime of the parent 
massive particle In particular, Holtmann and the present authors || performed the 

Maximum Likelihood analysis including both theoretical and observational errors and ob- 
tained the precise constraints. 

After Ref. || was published, several new observational data of light elements were re- 
ported. As for the 4 He abundance, it was still unclear whether the observational value of 
the primordial 4 He mass fraction Y is low (~ 0.234) |4],U or high (~ 0.244) ||. However, 
Fields and Olive considered the Hel absorption effect and reanalyzed the data , and they 
obtained a relatively middle value of Y (~ 0.238). On the other hand, as for the primordial 
D/H, although low values of D/H (~ 10~ 5 ) |§ had been measured and regarded as the 
primordial abundance, a relatively high value of D/H (~ 10 -4 ) was claimed again by Tytler 
et al. in the high redshift QSO absorption systems ||. In their paper they stressed that 
while the data may be inadequate to definitely conclude it to be precise value, there is still 
a possibility of the high D/H. 

On the theoretical side, recently it was claimed that the severest constraint on the 
radiatively-decaying particle may be from the non-thermal production of 6 Li, which is a 
secondary 6 Li production due to the background 4 He and the energetic T or 3 He produced 
by the 4 He photodissociationj] However, the observational data of the primordial compo- 
nent of 6 Li has large uncertainties. In addition, precise experimental data for the nuclear 
cross sections are not available. Therefore, it is unclear how important the non-thermal 6 Li 
production is once we take account of these uncertainties. 

With these new developments in theory and observation, we revise the previous con- 
straint on the radiative decay of long-lived particles. We obtain the photon spectrum by 
solving the Boltzmann equation numerically 0. In addition we perform the Monte Carlo 
simulation which includes both the experimental and theoretical errors. Then, we estimate 
the confidence levels by performing the Maximum Likelihood method including both the 



Such a possibility of the secondary process had already been pointed out by the earlier works for 
hadronic decaying particles 
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theoretical and the observational errors. 

This paper is organized as follows. In Sec. II we briefly review the current status of the 
observational data. In Sec. Ill we introduce the formulations for the photodissociation and 
non-thermal 6 Li production. In Sec. IV we compare the theoretical predictions with the 
observations. Sec V is devoted to conclusions. 



II. OBSERVATIONAL LIGHT ELEMENT ABUNDANCES 

Here we summarize the current status of the observational light element abundances. The 
primordial D/H is measured in the high redshift QSO absorption systems. Recently a new 



deuterium data was obtained from observation of QSO HS 0105+1619 at z = 2.536 fLzf . It 
was found that the cloud is neutral and has a simple structure. Five Lyman series transitions 
caused by D and H were observed there. The reported value of the deuterium abundance 
was relatively low, (D/H) obs = (2.54 ± 0.23) x 10~ 5 . Combined with the previous "Low D" 
data which were obtained by the clouds at z = 3.572 towards Q1937-1009 and at z = 2.504 
towards Q1009+2956 ||, the primordial abundance is obtained as 

LowD : (D/H) obs = (3.0 ± 0.4) x 10" 5 . (1) 

We call this value "Low D." On the other hand, Webb et al. oobserved high deuterium abun- 
dance in relatively low redshift absorption systems at z = 0.701 towards Q1718+4807 |p"3f . 

HighD : (D/H) o6s = (2.0 ± 0.5) x 10" 4 . (2) 

Tytler et al. @ also observed the clouds independently and obtained the similar value. Since 
Webb et al. and Tytler et al. did not obtain the full spectra of the Lyman series in their 
observations, the precise fitting of D /H based on the "High D" data might be inadequate. 
However, the possibility of "High D" have not been excluded yet. Therefore, we also consider 
the possibility of "High D" and include it in our analysis. 

For 3 He, we use the pre-solar measurements. In this paper, we do not rely upon models of 
galactic and stellar chemical evolution because of the large uncertainty in extrapolating back 
to the primordial abundance. But it is reasonable to assume that 3 He/D is an increasing 
function of the cosmic time, because D is the most fragile isotope and is always destroyed 
whenever 3 He is destroyed. Using the solar-system data reanalyzed by Geiss |14j, 



rfl & = ( 3 He/D)° &S = 0.591 ± 0.536, (3) 



where denotes the pre-solar abundance. We take this to be an upper bound on the 
primordial 3 He to D ratio r^f: 

r°s b 2 < (4) 

Although in the standard scenario the theoretical prediction satisfies the above constraint, 
4 He photodissociation produces both D and 3 He and can raise the 3 He to D ratio [15]. Hence, 
we include this constraint into our analysis. 

The primordial 4 He mass fraction Y is inferred from observation of recombination lines 
from the low metallicity extragalactic HII regions. Since 4 He is produced with Oxygen in 
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stars, the primordial value is obtained to regress to the zero metallicity O/H — > for the 
observational data. Recently, Fields and Olive [|7| reanalyzed the data including the Hel 
absorption effect and they obtained 



Y obs = 0.238 ± (0.002) stat ± (0.005)^, (5) 

where the first error is the statistical uncertainty and the second error is the systematic one. 
We adopt the above value as the observational Y. 

The primordial 7 Li/H is observed in the Pop II old halo stars. We adopt the recent 



measurements by Bonifacio and Molaro [jrj 

, obs 



log 
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( 7 Li/H)° S = -9.76 ± (0.012) stat ± (0.05) syst ± (0.3) add . (6) 



Here we have added the additional uncertainty for fear that 7 Li in halo stars might have 
been supplemented (by production in cosmic-ray interactions) or depleted (in stars) [|H]]. 

It is much more difficult to observe the primordial component of 6 Li because 6 Li is so 
much rarer than 7 Li. Unfortunately, enough data have not been obtained to find the "Spite 
plateau" of 6 Li. However, we can set an upper bound on 6 Li/ 7 Li, since it is generally believed 
that the evolution of 6 Li is dominated by the production through p,a-C,N,0 cosmic ray 
spallation (reactions of cosmic rays with the interstellar medium). Intrinsically the models 
of the nucleosynthesis through the cosmic ray spallation were motivated to simultaneously 
agree with whole the observational Li-Be-B abundances [|T8|-|2"0|| . On the other hand, recently 
it was claimed that the observational 6 Li abundance in halo stars is too abundant from the 
point of view of the cosmic ray energy if 9 Be is fit by the model of the cosmic-ray metal ||21|| . 
Therefore, there seems to be some uncertainties in the models of the cosmic ray spallation. 
In this situation, however, at least it would be safe to assume that 6 Li abundance increases 
as the metallicity increases. Today we observe only the 6 Li to 7 Li ratio in low-metallicity 
([Fe/H] < -2.0) halo stars H, 



r 



obs 

6,7, halo 



Ti/ 7 Li)°l = 0.05 ± 0.02. (7) 



We take this value as an upper bound on the primordial value r, 



obs 
6,7' 



r, 



" hs < rtlkaio- (8) 



6,7 — ' 6,7,/iaZo- 



III. PHOTODISSOCIATION AND NON-THERMAL PRODUCTION OF 6 LI 



A. Photodissociation 



In order to discuss the effect of high-energy photons on BBN, we need the shape of the 
photon spectrum induced by the primary high-energy photons from the decay of the massive 
particle X. In the thermal bath (mixture of photons 7bg, electrons e^ G , and nucleons Nbg), 
high energy photons lose their energy by the following cascade processes: 

• Double-photon pair creation (7 + 7bg — > e + + e~) 
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• Photon-photon scattering (7 + 7bg ^7 + 7) 

• Pair creation in nuclei (7 + Abg ~~ * e+ + e_ + N) 

• Compton scattering (7 + eg G — > 7 + e~) 

• Inverse Compton scattering (e ± + 7bg — > e ± + 7) 

In this study we numerically solved the Boltzmann equation including the above processes, 
and obtained the distribution function of photons, f y (E y ). 

The cascade photons induce the photodissociation of the light elements, which modifies 
the result of standard BBN (SBBN). The evolutions of the light nuclei abundances are 
governed by the following Boltzmann equation: 



dn N 

— — + 3Hn N 
at 



dn 



N 



J SBBN N 



dt 

+ J2 n N" / dE 7 a N u^ N (E 7 )fy(E 7 ), (9) 

AT"// M 



where n N is the number density of the nuclei N, and [drajv/^lsBBN denotes the SBBN 
contribution to the Boltzmann equation. In Table |, we list the photodissociation processes 
included in our computation. In this study the model parameters are the baryon to photon 
ratio (77), the lifetime of X (tx), the mass of X (mx), and the yield variable Yx of X after 
electron-positron annihilation, 

Y x = n x /n 7 , (10) 

where n 7 is the number density of photon.^ In this paper we assume that X decays only 
into photons, i.e., mxYx corresponds to Ap 7 /n 7 . Then, the photodissociation rates depend 
on the combination mxYx which characterizes the amount of the energy of the injected 



photons Ap 7 as far as mx is much larger than 20 MeV f31 |. 



B. Non-thermal 6 Li production 



As pointed out by Jedamzik [JTUJ , both T and 3 He are produced through the photodisso- 
ciation of 4 He, 

4tt , fn + 3 He ( , 

He + 7 ^| p + T (11) 

They are still energetic and have kinetic energies enough to produce 6 Li through the following 
processes with the background 4 He: 



Note that in the reference ||, Yx = nx/n^ is defined before electron-positron annihilation (e + e 



ann.). Then they have a relationship | after e+e-ann. = Tf Y x\beiore e+e-ann 
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T + 4 He — > 6 U + n, (12) 
3 He + 4 He — ► 6 Li + p, (13) 

until they are stopped by the ionization loss through the plasma excitation in the electro- 
magnetic plasma. The threshold energy of the 6 Li production is E$ Ue = 4.03MeV for 3 He, 
and = 4.80MeV for T. Then, the abundance of 6 Li produced through the non-thermal 
(NT) process in Eq. (|TJ) is governed by 



Li 



dt 



NT 



f°° f{E-y-El h )/A ( dE\ 

n4Re I h th ^7 cr4 He( 7lP )T (S 7 )/ 7 (£? 7 ) / n4 He CT T ( ain )6 Li (E) — dE, (14) 

where n6 Li (n4 He ) denotes the number density of 6 Li( 4 He). (T4 He ( 7;P ) T (^ 7 ) is the cross section 
of the 4 He photodissociation, Ef 1 is the threshold energy of the photodissociation process, 
fj(Ej) is the photon spectrum which are obtained by solving the Boltzmann equation, 
c"T(a,n) 6 Li(-E) is the cross section of the process in Eq. (JT3). dE/dx denotes the rate of the 
ionization loss while the charged particle T is running a distance dx in the electromagnetic 
plasma. The rate of the ionization loss is expressed by [|32| 

dE Z 2 a 2 [ Am e p 2 \ 

^ ln b^J' (15) 

where uj 2 is plasma frequency (= 4irn e a/m e ), n e is the electron number density, m e is electron 
mass, Z is the charge, A ~ 0(1) is a constant and (3 is the velocity. The effect of the process 
(|13|) is described by replacing the suffix T with 3 He in Eq. (|I4"D . 

We include the above two processes of the non-thermal 6 Li production in BBN code and 
compute the 6 Li abundance. In the computation we adopt the experimental cross section 
o"T(a,n) 6 Li = 35 ± 1.4 nib ^3[ commonly for the two processes. Because we have only one data 



point at the kinetic energy Et = 28 MeV in the laboratory system, we assume that the cross 
section is constant for whole the energy region and neglect its energy dependence. Then, 
we integrate the second factor in Eq. ( |I4"D up to a high energy. One can easily find that 
there exists a serious problem in this procedure if it is compared to the case of the original 
photodissociation where the photodissociation rates steeply decrease as the energy increases. 
Because we have the experimental data for the 4 He photodissociation rates only up to about 
100 MeV for the photon energy ]2"T|-f2"^], we should interpolate the photodissociation rates in 
a high energy region because of the mild dumping of the integrand. Then, the integration 
has a large uncertainty (~ 20 %) when we change the upper limit of the integration from 
500 MeV to 1 GeV.0 Therefore, in this situation we adopt 20% errors for the non-thermal 
6 Li production rates and perform the Monte Carlo simulation including them.0 



3 In addition, there may be another larger uncertainty which comes from the differences of the 
method for the interpolation because we do not know the correct shape of the cross sections. In 
this case, the obtained constraint would be weaker. 

4 If the cross section o"T(a,n) 6 Li decreases at high energy like other nuclear interactions, the 6 Li 



6 



C. Constraint from cosmic microwave background 

In addition to the photodissociation process, there also exists an another constraint. 
A radiative decay process releases a net photon energy into the electromagnetic plasma. 
The emitted photons should be thermalized soon, otherwise the photon spectrum deviates 
from the blackbody, which contradicts the observation of the cosmic microwave background 
(CMB) [p4|l - This leads to the following constraints: 

m x Y x < 2.0 x l(T 12 GeV (j^-) * , (16) 

\10 iU sec/ 

for /i-distortion (1.8 x 10 6 sec {Q B h 2 /0.02)^ < r x < 2.3 x 10 9 sec {Q B h 2 /0.02)), and 

m x Y x < 1.9 x 10- 12 GeV (j^-) * , (17) 

\10 iU sec/ 

for y-distortion (2.3 x 10 9 sec (tt B h 2 /0.02) < r x < 10 12 sec). 

IV. COMPARISON WITH OBSERVATIONAL LIGHT ELEMENT ABUNDANCES 

In Fig. [I] we plot the theoretically predicted 6 Li to 7 Li ratio (= r^) in (t x , m x Y x ) plane. 
The solid line represents the model parameters which predict the observational mean value 
of r l Q 7 and the dashed line denotes the observational 2-er upper bound. From the figure, one 
may think that the mean value of the theoretical prediction constrains m x Y x severely. We 
should bear in mind, however, that the theoretical prediction has a large uncertainty which 
comes from the errors of the production rates, and in addition the observational constraint 
also has a large error. To take account of these uncertainties systematically, we performed 
the Maximum Likelihood analysis P including both the theoretical and the observational 
errors. Here we assume that the theoretical predictions of (D/H) th , Y th , logio[( 7 Li/H)* /l ], 
r 32 = ( 3 He/D)* h , and 7g 7 obey the Gaussian probability distribution functions (p.d.f.'s) 
with the widths given by the l-er errors. Concerning the observational values, (D/H) ofes , 
Y obs , and logio[( 7 Li/H) ofes ] are assumed to obey the Gaussian p.d.f.'s while we treat r^l and 
rg & 7 as non-Gaussian variables 0. 

In Fig. |] we plot the results of the x 2 fitting by using the method of the Maximum 
Likelihood analysis. The solid (dashed) line denotes the Low D (High D) constraint. The 
dotted line denotes the upper bound from the CMB constraint. In the figure, the region 
below the lines are consistent with the observations. The constraint from the CMB is almost 
always weaker than that from BBN. The main feature of the difference between High D and 
Low D is that the Low D constraint is severer than High D for a relatively long lifetime case 
(tx ~ 3 x 10 6 sec). That is because High D constraint modestly allows the overproduction 
of 3 He accompanying the 4 He photodissociation. On the other hand, the High D constraint 



production is less important. As shown later, the resultant constraint is not changed even if we 
neglect the 6 Li production. 
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is more stringent for shorter lifetimes since the D dissociation is more important than the 
4 He photodissociation. 

The obtained upper bound does not change our earlier results so much ||. It became 
slightly weaker because we included the r\ dependence for the photodissociation rates (T oc 
I/77) in this analysis. f\ We find that the non-thermally produced 6 Li mildly contributes to 
the bound.0 The main reason is that both the theoretical computation and observational 
data have very large uncertainties which amount to about 30 - 40 %. 

Assuming that the parent massive particle is the gravitino and that it dominantly decays 
into a photon and a photino (^3/2 — > 7 + 7), the lifetime r 3 / 2 is related to the gravitino mass 
m 3/2 as 

r 3/2 ~ 4 x 10 5 sec x (m 3/2 /l TeV)" 3 . (18) 

Assuming that the gravitino is produced through the thermal scattering in the reheating 
process after inflationfj we relate the abundance Y 3 / 2 = n 3 / 2 /n 7 of the gravitino with the 
reheating temperature T R 0, 

V 3/2 ~ 1.1 x 10~ n x (2Vl0 10 GeV). (19) 

In Fig. § we plot the upper bound on the reheating temperature after inflation at 95% 
C.L. as a function of the gravitino mass. Here we can read off the constraint by using 
the relationship of the scaling, Ap 7 /n 7 = |m. 3 / 2 Y" 3 / 2 (= m x Y x ) because we assumed X 
decays into two photons. From the figure we can obtain the upper bound on the reheating 
temperature: 



m 3/2 = 


100 GeV 


( r 3/2 


~ 4 x 


10 8 


sec) 


:T R 


< 


1 x 10 7 GeV, 


m 3 / 2 


= 1 TeV 


( r 3/2 


~ 4 x 


10 5 


sec) 


:T R 


< 


1 x 10 9 GeV, 


m 3 / 2 


= 3 TeV 


( r 3/2 


~ 1 X 


10 4 


sec) 




< 


9 x 10 11 GeV, 



at 95% C.L. 



5 The rj dependence is understood as follows. The soft photons produced in the electromagnetic 
cascade scatter off the background electrons and nucleons and lose their energy. Thus the number 
density of soft photons with energy larger than the threshold decreases as scattering rate which is 
proportional to r]. Therefore, the photodissociation rates are proportional to 1/rj. 

6 Tritium is unstable with the lifetime tt = 5.614 x 10 8 sec, and decays into 3 He whose charge is 
two. Thus, because 3 He subjects to stop easier than T by the ionization loss, we might overestimate 
the 6 Li production in parameter regions where the stopping time T s ^ p — $z(dE/dt)- 1 dE ~ 2.5 x 
10 9 sec(T/eV) _3 (.E/MeV) 3//2 is longer than the lifetime of Tritium, i.e., for T < 1.7eV. Therefore, 
at a long lifetime tx ^ 5 x 10 11 sec, our constraint might become weaker by about factor two. 
However, it is expected that the effect would not change the result significantly because the 3 He 
overproduction gives a severer constraint there. 

"Although these days it was claimed that gravitinos are also produced in the preheating epoch 



non-thermally p5|-37|, we do not consider such processes here because there are some ambiguities 
on the estimations and they depend on various model parameters. If the non-thermal production 
is effective, however, the obtained constraint might be severer. 
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V. CONCLUSION 



In this paper we have studied the effects on primordial nucleosynthesis of the radiative 
decay of a long-lived massive particle X using new observational data. We have also con- 
sidered the non-thermal 6 Li production caused by energetic T and 3 He produced by the 4 He 
photodissociation. We obtained the photon spectrum through the electromagnetic cascade 
process by solving a set of Boltzmann equations numerically. In addition, to estimate the 
theoretical errors we performed Monte Carlo simulation including the theoretical uncertain- 
ties which come from those of nuclear reaction rates. To obtain the degree of agreements 
between theory and observation, we performed the Maximum Likelihood method and the 
X 2 fitting including both the observational and theoretical errors. 

As a result we have obtained the upper bound on the abundance mxYx as a function of 
its lifetime tx- The result does not change our previous works significantly. This is because 
the theoretical and observational errors for 6 Li are significantly large, and it contributes to 
the constraints more weakly than the 3 He overproduction accompanying the 4 He photodis- 
sociation. Therefore, we have concluded that it is premature to emphasize the importance 
of the non-thermal production of 6 Li. 

We have also applied the results obtained by a generic radiatively decaying particle 
to gravitino ^3/2, an d we have got the upper bound on the reheating temperature after 
primordial inflation as a function of the mass, T R < 10 7 — 10 9 GeV for m 3 / 2 = 100 GeV — 
1 TeV (95 % C.L.). 

VI. ACKNOWLEDGMENTS 

The work of M.K. is supported by Priority Area 707 "Supersymmetry and Unified The- 
ory of Elementary Particles." The work of T.M. is supported by the Grant-in-Aid for Sci- 
entific Research from the Ministry of Education, Science, Sports and Culture of Japan No. 
12047201. 



9 



REFERENCES 



[1] J. Ellis, G.B. Gelmini, J.L. Lopez, D.V. Nanopoulos, and S. Sarkar, Nucl. Phys. B373, 
399 (1992). 

[2] M. Kawasaki and T. Moroi, Prog. Theor. Phys. 93 879 (1995). 

[3] E. Holtmann, M. Kawasaki, K. Kohri, and T. Moroi, Phys. Rev. D60 023506 (1999). 

[4] B.E.J. Pagel et al, Mon. Not. R. Astron. Soc. 255 (1992) 325. 

[5] K.A. Olive, G. Steigman, and E.D. Skillman, Astrophys. J. 483 (1997) 788. 

[6] Y.I. Izotov, T.X. Thuan, and V.A. Lipovetsky, Astrophys. J. Suppl. Series, 108 (1997) 

1; Y.I. Izotov and T.X. Thuan Astrophys. J., 500 (1998) 188. 
[7] B.D. Fields and K.A. Olive, Astrophys. J., 506 177 (1998). 
[8] S. Buries and D. Tytler, Astophys. J. 507 732 (1998). 

[9] D. Tytler, S. Buries, L. Lu, X-M. Fan and A. Wolfe, Astron. J. 117 (1999) 63. 
[10] K. Jedamzik, Phys. Rev. Lett. 84 (2000) 3248. 

[11] S. Dimopoulos, R. Esmailzadeh, L.J. Hall, and G.D. Starkman, Astrophys. J. 330 545 

(1988); Nucl. Phys. B311 699 (1989). 
[12] J.M. O'Meara, D. Tytler, D. Kirkman, N. Suzuki, J.X. Prochaska, D. Lubin, and A.M. 

Wolfe, |astro-ph/00lTT79] 
[13] J. K. Webb et al, Nature 388 (1997) 250. 

[14] J. Geiss, in Origin and Evolution of the Elements, edited by N. Prantzos, E. Vangioni- 
Flam, and M. Casse (Cambridge University Press, Cambridge, 1993) 89. 

[15] G. Sigl, K. Jedamzik, D.N. Schramm, and V.S. Berezinsky, Phys. Rev. D52 6682 
(1995), |aito-ph/ 9503094 . 



[16] P. Bonifacio and P. Molaro, Mon. Not. R. Astron. Soc. 285 847 (1997). 
[17] B.D. Fields, K. Kainulainen, K.A. Olive, and D. Thomas, New Astron. 1 (1996) 77, 
^stro-ph/96030m 



[18] D.K. Duncan, D.L. Lambert, and M. Lemke, Astrophys. J. 401 584 (1992); M. Casse, 
R. Lehoucq, and E. Vangioni-Flam, Nature 373 318 (1995); R. Ramaty, B. Kozlovsky, 
and R.E. Lingenfelter, Astrophys. J. Lett.438 L21 (1995). 

[19] M. Lemoine, D.N. Schramm, J.W. Truran, and C.J. Copi, Astrophys. J. 478 554 (1997). 

[20] B.D. Fields and K.A. Olive, New Astron. 4 255 (1999). 

[21] R. Ramaty, S.T. Scully, R.E. Ligenfelter, and B. Kozlovsky, Astrophys. J. 534 747 
(2000). 

[22] V.V. Smith, D.L. Lambert, and PE. Nissen, Astrophys. J. 408 262 (1993); L.M. 
Hobbs and J. A. Thorburn, Astrophys. J. 491 772 (1997); V.V. Smith, D.L. Lambert, 
and P.E. Nissen Astrophys. J. 506 923 (1998); R. Cayrel et al, Astron. Astrophys. 343 
923 (1999). 

[23] R.D. Evans, The Atomic Nucleus (McGraw-Hill, New York, 1955). 
[24] R. Pfeiffer, Z. Phys. 208 (1968) 129. 

[25] D.D. Faul, B.L. Berman, P. Mayer, and D.L. Olson, Phys. Rev. Lett. 44 (1980) 129. 

[26] A.N. Gorbunov and A.T. Varfolomeev, Phys. Lett. 11 (1964) 137. 

[27] J.D. Irish, R.G. Johnson, B.L. Berman, B.J. Thomas, K.G. McNeill, and J.W. Jury, 

Can. J. Phys. 53 (1975) 802. 
[28] C.K. Malcom, D.V. Webb, Y.M. Shin, and D.M. Skopik, Phys. Lett. B47 (1973) 433. 
Yu.M. Arkatov, P.I. Vatset, V.I. Voloshchuk, V.A. Zolenko, I.M. Prokhorets, and V.I. 
ChimiP, Sov. J. Nucl. Phys. 19 (1974) 598. 



10 



[30] B.L. Berman, Atomic Data and Nuclear Data Tables 15 (1975) 319. 

[31] M. Kawasaki and T. Moroi, Astrophys. J. 452, 506 (1995). 

[32] J.D. Jackson, Classical Electrodynamics (Wiley, New York). 

[33] J. A. Koepke and R.E. Brown, Phys. Rev. C16 (1977) 18. 

[34] Fixsen et al., Astrophys. J., 47 3 576 (1996). 

[35] R. Kallosh, L Kofman, A Linde, and A.V. Proeyen, Phys. Rev. D61 103503 (2000). 

[36] G.F. Giudice, A. Riotto, and I. Tkachev, JHEP 9908 009 (1999). 

[37] D.H. Lyth, Phys. Lett. B476 356 (2000). 



11 



TABLES 



Photodissociation Reactions 



1-d Uncertainty Threshold Energy 



Ref. 



1. 


D + 7 ^ 


p + n 


6% 


2.2 MeV 


2. 


T + 7 ^ 


n + D 


14% 


6.3 MeV 


3. 


T + 7 ^ 


p + 2n 


7% 


8.5 MeV 


4. 


3 He + 7 


^p + B 


10% 


5.5 MeV 


5. 


3 He + 7 


— > n + 2p 


15% 


7.7 MeV 


6. 


4 He + 7 


^p + T 


4% 


19.8 MeV 


7. 


4 He + 7 


—> n + 3 He 


5% 


20.6 MeV 


8. 


4 He + 7 


— > p + n + D 


14% 


26.1 MeV 


9. 


6 Li + 7 - 


-> anything 


4% 


5.7 MeV 


10. 


7 Li + 7 - 


-> 2n + anything 


9% 


10.9 MeV 


11. 


7 Li + 7 - 


-> n + 6 Li 


4% 


7.2 MeV 


12. 


7 Li + 7 - 


-> 4 He + anything 


9% 


2.5 MeV 


13. 


7 Be + 7 


6 Li 


4% 




14. 


7 Be + 7 


— > anything except 6 Li 


9% 





[23] 



m 
m 



1 27||28 1 



m 
H3on 



TABLE I. List of photodissociation processes, and the 1-a uncertainty in the cross sections. 
Since there is no experimental data on photodissociation of 7 Be, we assume that the rate, threshold, 
and uncertainty for Reaction 13 is the same as for Reaction 11, and the rate for Reaction 14 is the 
sum of the rates for Reactions 10 and 12. 
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FIGURES 
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FIG. 1. Plot of 6 Li to 7 Li ratio in (tx, mxYx) plane for various baryon to photon ratio 
(77 = n B /n 7 ) in (a) rj = 2 x 1CT 10 , (b) rj = 4 X 1CT 10 , (c)i) = 5x 1CT 10 , and (d) 77 = 6 x 1CT 10 . 
The solid line denotes the observational mean value of 6 Li / 7 Li and the dashed line denotes the 
observational 2-a upper bound. 
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t x (sec) 



FIG. 2. Plot of the contour of the confidence level in {T x ,mxYx) plane. The solid (dashed) 
line denotes the 95% C.L. for Low D (High D) projected on rj axis. The dotted line denotes the 
upper bound which comes from CMB constraint. 
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FIG. 3. Plot of the contour of the confidence level in (7713/2, Tr) plane. The solid (dashed) line 
denotes the 95% C.L. for Low D (High D). The dotted line denotes the upper bound which comes 
from CMB constraint. 
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